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Abstract: We consider the ratio of cross sections of double-to-single Higgs boson 
production at the Large Hadron Collider at 14 TeV. Since both processes possess 
similar higher-order corrections, leading to a cancellation of uncertainties in the ra- 
tio, this observable is well-suited to constrain the trilinear Higgs boson self-coupling. 
We consider the scale variation, parton density function uncertainties and conserva- 
tive estimates of experimental uncertainties, applied to the viable decay channels, to 
construct expected exclusion regions. We show that the trilinear self-coupling can 
be constrained to be positive with a 600 fb^^ LHC dataset at 95% confidence. More- 
over, wc demonstrate that we expect to obtain a +30% and —20% uncertainty on 
the self-coupling at 3000 fb^^ without statistical fitting of differential distributions. 
The present article outlines the most precise method of determination of the Higgs 
trilinear coupling to date. 
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1. Introduction 

One of the aims of the Large Hadron Collider (LHC) is to search for the agent of 
electroweak symmetry breaking (EWSB), which in its minimal form is the Standard 
Model (SM) Higgs boson (H). Recently, both the ATLAS and the CMS collabora- 
tions have observed a new state with a mass of about 125 GeV, whose properties are 
in substantial agreement with the SM Higgs boson [1-5]. The quest for understand- 
ing the mechanism behind EWSB docs not end with the discovery of this particle. 
It is crucial to test the Higgs sector to its full extent, measuring the couplings of 
the Higgs boson to gauge bosons and matter fields [6-30], and also to probe its 
self-interactions [31-35]. After EWSB, the Higgs potential can be written as 



In the SM, Ag^^ = X^/^hh = {Mh/2v^) ~ 0.13 for a Higgs mass of Mh ^ 125 GeV 



and a vacuum expectation value of f ^ 246 GcV. Wc can also define normalized 
couplings A = Xhhh / X^HHH ^^^d A = Xhhhh / Xf^HHH- With an extended Higgs 
sector, as is common in many new physics models beyond the SM, these couplings 
will deviate from the SM values. Therefore, measuring these two couplings is crucial 
to reveal the true nature of the Higgs boson. At the LHC, the quartic coupling A 
may be probed via triple Higgs boson production. However, its tiny cross section [36] 
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makes it very difficult, if not impossible, to do so. On the other hand, the trilinear 
coupling A can be measured in Higgs boson pair production, pp HH, which may 
be discovered at a large luminosity phase of the LHC. 

The discovery potential for Higgs boson pair production at the LHC has been 
studied in [32-35,37]. In Refs. [32,37], constraints were placed on A using statistical 
fits to the shape of the visible mass distributions of the final decay products of the 
Higgs pairs, whereas Refs. [33,34] focused on the estabhshment of the Higgs pair pro- 
duction process using improved techniques originating mainly from developments in 
the understanding of boosted jet substructure [38,39]. In Ref. [35] the final state 6677 
was revisited as well as bbr^r^ and bbW^W^ (fully leptonic), without making use 
of jet substructure techniques (although boosted Higgs bosons were required). The 
present article concentrates on using the results from the available phenomenological 
studies along with the best available theoretical cross section calculations and con- 
servative estimates of the experimental uncertainties, to demonstrate the possibility 
of constraining the trilinear Higgs self-coupling at the LHC. 

The article is organised in the following way: in Section 2 we dissect the Higgs 
boson production cross sections and in Section 3 we examine the theoretical uncer- 
tainties on the ratio of cross sections of double-to-single Higgs production. Then, in 
Section 4, we present the expected constraints obtained at integrated luminosities of 
600 fb~^ and 3000 fb~^. We conclude in Section 5. 

2. Dissection of the cross sections 

The Higgs boson pair production cross section is dominated by gluon fusion, as is 
the single production cross section [40,41]. For the pair production, other modes, 
like qq — )■ qqHH,VHH, tiHH are a factor of 10-30 smaller [42,43], and thus we 
do not consider them in the rest of our analysis. At leading order (LO), there are 
two main contributions: a diagram containing a 'triangle' loop, and one containing 
a 'box' loop of top (or bottom) quarks, as shown in Fig. 1. The production of a 
single, on-shell Higgs boson only contains a diagram of the 'triangle' type. The 
triangle diagram can only contain initial-state gluons in a spin-0 state, whereas the 
box contribution can contain both spin-0 and spin-2 configurations. Therefore, there 
are two Lorentz structures involved in the box diagram matrix element. At LO, we 
may write, schematically: 

'y'H'l^-WrC^}+PrC^:h' + ll\C':L\\ (2.1) 

where C^} represents the matrix element for the triangle contribution and C^l-^ 
represents the matrix element for the two Lorentz structures (i = 1,2) originating 
from the box contribution [40,44]. The parameters ai, ^1 and 71 are given in terms 



- 2 - 



of the Standard Model Lagrangian parameters, by: 



A = 71 = Vt , (2.2) 

where A is the (normahzed) Higgs triple coupling defined in the previous section 
and yt is the Standard Model Hti coupling (as defined after electroweak symmetry 
breaking and assumed to be real) normalized to the SM value. ^ In contrast, the 
single Higgs cross section, again, schematically, will only contain the matrix element 
squared IC^^^p- 




H 



Figure 1: The Higgs pair production diagrams contributing to the gluon fusion process 
at LO are shown. 

We have performed numerical fits using the results of the hpair program [45], 
used to calculate the total cross section for Higgs boson pair production at lead- 
ing and approximate next-to-leading (NLO) orders. The fits were done employing 
MSTW20081o68cl and MSTW2008nlo68cl parton density functions [46] and using 
top and bottom quark masses of 174.0 GeV and 4.5 GeV respectively. We have 
obtained: 

a^o = 5.15A2y2 _ 25.1Xy^ + 38.1yt , 

a^^o = 9MX'y^ - 46.9Ayf + 71.6yt . (2.3) 

It is evident from Eqs. (2.1) and (2.3) that the Higgs pair production cross 
section contains an interference term proportional to {Xyf)- Hence, for positive values 
of (Xyf) the cross section is reduced, whereas for negative values, it is enhanced. The 
box squared term is dominant, and scales as yf, whereas the triangle squared term is 
subdominant due to the off-shell Higgs boson which then decays to the Higgs boson 

^In the previous discussion, Eqs. (2.1) and (2.2), assume implicitly that the contribution from 
bottom loops are neglible. In the SM case they contribute to the total cross section less than 
0.2%. We have checked numerically that the variation of the bottom Yukawa coupling effect on 
each form factor is less than 1%. In the rest of this article we include the bottom loops in our 
numerical calculations assuming that the bottom Yukawa coupling has the SM value. For the sake 
of simplicitly we assume the validity of Eqs. (2.1) and (2.2) in what follows. 
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pairs, and scales as A^y^. Also note that there exists a minimum value of cr^^ at 
A ~ 2A6yt. 

We note that the structure and hence the different contributions to the cross 
section can of course be modified if new physics that allows new resonances to run 
in the triangle and box loops is present [47-51]. Here, we are assuming that A and 
lit can be modified, but no other particles enter the diagrams. 

3. Ratios of cross sections 

It has been pointed out in Ref. [26] that the ratio of cross sections between Higgs 
pair production and single Higgs production: 



could be more accurately determined theoretically than the Higgs-pair production 
cross section itself. This is based on the fact that the processes are both gluon- 
initiated and the respective higher-order QCD corrections could be very similar. 
Hence, it is assumed that a large component of the QCD uncertainties drop out in 
the ratio Chh- Moreover, experimental systematic uncertainties that affect both 
cross sections may cancel out by taking the ratio. An example is the luminosity 
uncertainty, which should cancel out provided the same amount of data is used in 
both measurements. 

Here we investigate the extent to which the above assumptions are correct, using 
the available calculations for the cross sections. We start by considering the LO and 
NLO calculations for a{gg — HH) and a{gg — ?> H) at the LHC at 14 TeV. Using 
the MSTW20081o68cl and MSTW2008nlo68cl parton density functions [46], we show 
in Figs. 2 and 3 the cross sections as well as their ratios, Chh, as a function of 
the Higgs mass at both LO and NLO.^ We present the scale uncertainty obtained 
by varying the factorisation and renormalization scales (set to be equal) between 
[0.5/xo) 2.0/io]; where //q = Mh for the higlu program, used to obtain the single 
Higgs cross sections [53] , and //q = Mhh for the hpair program, used for the Higgs 
pair production cross sections [45]. The scale choices are the natural ones for each 
of the processes but we verified that the conclusions are not altered substantially by 
changing the hpair scale, i.e. the numerator, to equal the scale that appears in the 
denominator, yUg = Mh- Implicit in the calculation of the scale uncertainty of the 
ratio Chh, is the fact that the scale variation of the single and double Higgs cross 
sections between 0.5/io and 2.0yLio is fully correlated: i.e., we obtain the upper and 
lower variations of the ratio by dividing the cross sections with the same magnitude 

^It is important to note that the NLO calculation has been performed in the heavy top mass 
limit, and hence in the case of Higgs boson pair production, it is expected to be approximate. At 
LO, the accuracy of the large top mass approximation is 0(10%) [31,50,52]. 



a{gg HH) ^ crgjj 
(^{99 H) ~ (Th 
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of variation of the scale. This is an approximation that is justified since the two 
processes possess similar topologies, and is in fact one of the main insights in favour 
of using Chh- We also show, in the ratio, the resulting PDF uncertainty, calculated 
using the MSTW2008nlo68cl error sets according to the prescription found in [54]. 



/^ = 14TeV, 0.5/Jn </J<2.0/Jn, LP (MSTW2008lo68cl) 




120 121 122 123 124 125 126 127 128 

(GeV) 



Figure 2: The cross sections for single and double Higgs boson production at leading 
order using the MSTW20081o68cl PDF set. In the lower plot, the fractional uncertainty 
due to scale variation is shown in the blue band, as well as the PDF uncertainty in the 
green band. 

Several observations on the behaviour of the Chh ratio can be made. First of 
all, it is evident that the fractional uncertainty due to scale variation is reduced with 
respect to the individual calculations in both leading and next-to-leading orders: 
for the LO case, the individual cross sections have a ~ ±16% (single Higgs boson 
production) and ~ ±25% (double Higgs boson production) scale uncertainty, whereas 
the ratio has a ~ ±9% scale uncertainty. For the NLO case, it is reduced from 
~ ±17% (single and double Higgs boson production) to ~ ±1.5% for the ratio. ^ 

Furthermore, we can explicitly see that the uncertainty due to the QCD correc- 
tions partially cancels out: even though the individual K-f actors in the cross sections 
aH and auH are large, they are also very similar, both being ~ 2. As a consequence, 
the central value of the ratio only decreases by a small amount from ~1.25 to ~1.0. 
This is an indication that higher order corrections are quite likely to change the ratio 
by an even smaller fraction than the change from LO to NLO, when it is consid- 
ered at NNLO, whereas the single cross section has a K-factor of about ~1.5 [59] 

■^Note that in Ref. [55], threshold resummation effects in SM Higgs pair-production in soft- 
coUinear effective theory were considered. The authors claim a reduction of the scale uncertainty 
to 3%. For other resummation studies in single Higgs production see, for example [56-58]. 
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v/T = 14TeV, 0.5/Jn <n<2.0iin, NLO (MSTW2008nlo68cl) 
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Figure 3: The cross sections for single and double Higgs boson production at next-to- 
leading order using the MSTW2008nlo68cl PDF set. In the lower plot, the fractional 
uncertainty due to scale variation is shown in the blue band, as well as the PDF uncertainty 
in the green band. 

with respect to the NLO calculation.^ These findings support the idea of employing 
the fully correlated scale variation described before as a realistic estimate for the 
theoretical error. ^ 

The PDF uncertainties for the cross sections themselves are not shown since they 
are of the order of a few % and hence subdominant. The PDF uncertainty is also 
sub-dominant in the case of the LO ratio, as shown in Fig. 2. In the case of the NLO 
ratio, the PDF uncertainty becomes comparable to the scale uncertainty as can be 
seen in Fig. 3. Combining the two errors in quadrature would induce an error of 
±0{3%), still smaller than the ~ ±17% error on the NLO Higgs pair production 
cross section. To remain conservative, we will assume that the theoretical errors on 
Chh and ann are ±5% and ±20%, respectively, in what follows. 

4. Constraining the self-coupling 

In the studies conducted in Refs. [32,37], the Higgs self-coupling was constrained 
using the final states 6677, bbfi'^fi~ and W^W~W^W~ (in the high Higgs mass 
region). The constraints were obtained by fitting the visible mass distributions in 
each process for the signal and backgrounds. 

Here we chose to follow a different strategy: taking into account the facts that 
the different signal channels possess a relatively low number of events and that the 

^An equivalent calculation at NNLO does not presently exist for Higgs pair production. 
^Note that a detailed study of the scale and parton density function uncertainties in Higgs pair 
production can be found in Rcf . [35] . 
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shapes of distributions for the backgrounds (and even the signal) are not always very 
well known, we employ only information originating from the rates. Furthermore, 
we use the theoretically more stable ratio Chh between the double and single Higgs 
production cross sections, examined in the previous section. We focus on luminosities 
of 600 fb~^ and 3000 fb~^ that can be respectively obtained by ATLAS and CMS 
together in the first long-term 14 TeV run, or by the individual experiments in an 
even longer-term run at the same energy. We do not attempt to combine between the 
individual channels, as this will require a more detailed study from the experimental 
collaborations. 



4.1 Variation with self-coupling and top quark Yukawa 

To quantify the possible region that can be constrained using the ratio Chh^ we first 
examine the behaviour of the cross section for Higgs pair production and the ratio 
Chh at 14 TeV, when varying the self-coupling A, as well as the top Yukawa, yt- It 
is important to consider variation of the top quark Yukawa determination, since the 
production rates of both double and single Higgs production can be substantially 
affected. Moreover, the expected accuracy on the top quark Yukawa is expected to 
be ±C(15%) at 300 fb^^ of LHC data at 14 TeV [60]. 

We show the cross section (Jhh and ratio Chh at = 1 as a function of A, as 
well as both quantities at A = 1 as a function of yt in Figs. 4 and 5, respectively. 
Evidently, the effects of both A and yt are significant: the cross section varies from 
~ 30 fb at (A,?/i) = (1, 1) (i.e. the SM values) to ~ 125 fb at [\yt) = (-1, 1) and 
~ 300 fb at {\,yt) = (1,1-6). The ratio itself varies from ~ 1 at {\,yt) = (1,1) to 
~ 3.5 at {\,yt) = (—1, 1) and {X,yt) = (1, 1-6). It is obvious that negative values of 
A can be excluded sooner than the positive values, since the cross section and ratio 
of cross sections both increase fast with decreasing A. 



S/f =1 Mg =125 GeV, LHC(B14 TeV, MSTW2008nlo 68cl =1 , Mg =125 GeV, LHC(ai4 TeV, MSTW2008nlo 68cl 
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Figure 4: The cross section for double Higgs production and the ratio Chh at next-to- 
leading order using the MSTW2008nlo68cl PDF set, as a function of A at = 1. 



We note that negative values of yt are currently viable [21] and physical, and 
could arise in beyond-the-SM physics models. Since Higgs pair production only 
depends on the sign of the product (Xyt), the corresponding values for yt < 0, \ > 



-7 - 



A=l Mjj =125 GeV, LHC(gil4 TeV, |VISTW2008nlo 68cl 

350 
300 
[250 
, 200 
' 150 



A = l 




Mg =125 GeV, LHC(g)14 TeV, |VISTW2008nlo 68cl 




Figure 5: The cross section for double Higgs production and the ratio Chh at next-to- 
leading order using the MSTW2008nlo68cl PDF set, as a function of yf at A = 1. 



are equivalent to those for the points with the same absolute values of the parameters 
but > 0, A < 0.6 

4.2 Deriving constraints 

The ratio Chh can be used to derive the expected constraints that can be obtained 
at a 14 TeV LHC Certain assumptions on the systematic uncertainties need to be 
made for the branching ratios related to each mode. We first define the following 
quantities: 



a^^lf = aHH X BR{hh) x BR(xx) 
a^^-anx BR{bb) , 



(4.1) 



where xx denotes the H ^ xx decay mode in question. Hence, we can invert the 
above relations to obtain: 



/^exp. 



a 



bbxx 
HH 



CT^ X BR{xx) 



(4.2) 



exp. 



which is the experimental measurement of the theoretical quantity Chh- 

Since the scope of this article is not a detailed experimental study, we now make 
several assumptions on the measurement uncertainties for each of the quantities in 
the ratio of Eq. (4.2). We focus on the region A G (—1.0,2.0), since the cross sec- 
tion becomes too small for A > 2 J According to Ref. [63], the branching ratio of 
H ^ hh times the cross section for single Higgs is expected to be known to ±20% 



^Note that the degeneracy that appears in Higgs pair production may be resolved through 
the study of different processes long before the Higgs self-coupling is probed. See, for example, 
Refs. [61,62]. 

^Note that there exists a minimum of the Higgs pair production cross section at A ^ 2.46 
(for yt — 1). Beyond this, the cross section values are symmetric with respect to values below 
A ~ 2.46. However, values of the coupling above this region can be considered to be 'large' from 
the perturbative point of view, and we do not consider them here. 
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after 300 fb~^ of data at 14 TcV, and hence we assume that the uncertainty on a^^ 
is ±20%. Similarly, according to [63], the uncertainties on BR(r"'"r~), BR(W^"'"W^~) 
and BR(77) are expected to be ±12%, ±12% and ±16%, respectively, at 300 fb~^. 
To remain conservative, we assume that going beyond 300 fb~^ of luminosity, there 
will be no improvement on these uncertainties. This can be true, for example, if the 
measurements are dominated by systematic uncertainties that cannot be improved 
further. Moreover, the uncertainty on the cross section of the measured final state, 
AcT^lf , is estimated by assuming that the Poisson distribution of the obtained num- 
ber of events can be approximated by a Gaussian, for simplicity. Hence, if we expect 
a number of B background events and we experimentally measure events, the 
error on the signal estimate, S = N — B, is given by AS* = \/N + B. The expected 
number of events for the studies we consider below were taken from [33,34,37]. We 
combine all the estimates of the uncertainties in quadrature for each mode to obtain 
an estimate of the total uncertainty: 

(AChhV _ (AafA' f ABR{xx) y ( Aa'^V 

\Chh J [ ) ^ \ BR(xx) ) ^ \ a^^ ) ' ' 

In what follows we also add the theoretical error estimates in quadrature to the 
above. 

As a first approximation we assume that the top quark Yukawa coupling has the 
SM value and is determined without error, i.e., that yt — 1. Thus, we consider a 

situation where the Standard Model is valid everywhere but in the Higgs potential. 
Inevitably, this introduces a degree of model-dependence on the derived A exclusion. 
We will consider relaxing the yt = I assumption below. Note that the possibility of 
yt = —1 can be studied simply by reflection: A —A. 
One can form the following question: 

Given an assumption for the 'true' value of the Higgs self-coupling, Atruej 
what is the constraint we expect to impose on A? 

The question is more elaborate than the one that was being asked in searches of the 
Higgs boson itself, as in that situation the absence of the Higgs Boson at a certain 
mass resulted in a specific prediction for the background rate. In the case of A, 
we have a continuous variable that has to be assessed, and the same experimental 
results could be due to different true values of A. The answer to the above question 
results in an 'exclusion' plot, calculated by drawing the curves that result in expected 
measurements that arc one or two standard deviations away from the central value 
which is assumed to be equal to that given by Atrue- By virtue of this definition, it 
is obvious that the central value itself is, of course, not expected to be excluded. 

Using Chh we draw such curves for 600 fb~^ of data in Figs. 6, 7 and 8 for the 
final states hhr^T^ , hhW'^W^ and 6677, respectively. To bring the three channels 
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to an equal footing, wc have rescaled the bbT^T~ cross section in [33] by employing 
a factor of 32.4/28.4 accounting for the central value of the NLO production cross 
section used in [34], and moreover, rescaled by 0.7^/0.8^ for a reduced r-jet tagging 
efficiency. For the bbW~^W~ mode in [34] we use the tauonic decays of the W 
bosons, and for the 6677 result in [37] we average between the 'hi' and 'lo' LHC 
results to get 6 versus 12.5 events at 600 fb"-*^. We have not rescaled the 6677 
analysis, since this was done for a Higgs of mass 120 GeV in [37]. In the lower 
panel of Fig. 6 we also show the exclusion regions extracted by using the Higgs 
pair production cross section measurement itself, with an associated uncertainty of 
±20%. Wc assume that the uncertainty on BR(66) is the same as that on o"^, namely 
±20%. It is obvious that the exclusion obtained from the cross section is expected to 
be weaker than that obtained by the ratio, due to the larger theoretical systematic 
uncertainty on the cross section itself. Moreover, the expected exclusion from ann 
will be more affected by experimental systematic uncertainties which would add to 
the errors. For completeness, we show the estimated fractional uncertainty on the 
ratio, A-Chh/Chh, used to extract the exclusion regions, for the different processes 
and investigated luminosities in Table 1. At high luminosity the uncertainties all tend 
to similar numbers since we have assumed that the other contributing uncertainties 
(ABR(a;a;) and Acr^) do not improve and they become systematic-dominated. These 
values are provided for completeness, as an indication, and merit further investigation 
by the experimental collaborations. 



Process 


^C„h/Chh (600 fb-i) 


^Chh/Chh (3000 fb-i) 


bbT^T~ 


0.335 


0.257 


bbW+W- 


0.558 


0.325 


bbjj 


0.738 


0.401 



Table 1: The fractional uncertainties on the ratio of double-to-single Higgs boson produc- 
tion cross sections, ^Chh/Chh, for the different channels and the two investigated LHC 
luminosities, 600 fb~^ and 3000 fb~^, using Mh = 125 GeV. These include the theoretical 
error due to the scale/parton density functions uncertainties, assumed to be 5%. 

The interpretation of the 'exclusion' curves is simple: as an example, if we assume 
or believe that the 'true' value of the triple Higgs coupling is Atrue = 1, then by 
examining Fig. 6 for the bbr^r^ mode at 600 fh^^, we can conclude that the expected 
experimental result should lie within A G (0.62, 1.52) with ~68% confidence (lo"), 
and A G (0.31, 3.08) at ~95% {2a) confidence. Wc expect to exclude any values 
outside this range after 600 fb~^, given the value Atrue = 1- We show the collected 
exclusion limits for Atrue — 1 ^-nd yt,tTue — 1 (i-e. the SM values) at la and 2a at 
600 fb"^ as well as the end-of-run LHC integrated luminosity of 3000 fb~^ in Table 2. 
The 3000 fb~^ values have also been calculated by assuming no improvement in the 
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uncertainty estimates that we have assumed at 600 fb^^. The table demonstrates an 
important conclusion: it is possible, using the discovery of the three viable channels, 
to constrain the trilinear coupling A to be positive at 95% confidence at 600 fb~^ 
(provided yt,true = !)• Moreover, a naive combination of the 'uncertainties', at la 
about Atrue, over the three channels indicates that a measurement of accuracy +30% 
and —20% is possible simply by using the rates at 3000 fb~^. 



Process 


600 fb-i {2a) 


600 R-)-^ (la) 


3000 fb-i 2a 


3000 fb-i la 




(0.31, 3.08) 


(0.62, 1.52) 


(0.45,1.99) 


(0.70, 1.37) 


bbW+W- 


(-0.04, 19.37) 


(0.41, 2.20) 


(0.33, 2.88) 


(0.63,1.50) 




(-0.29, 65.56) 


(0.26, 3.99) 


(0.20, 5.22) 


(0.56, 1.66) 



Table 2: The expected limits on A at la and 2a confidence levels, provided that Atrue and 
yt,true have their SM values: Atrue = 1, yt,true = 1- The results have been derived using 
Chh and are shown for 600 fb"^ and 3000 fb^^ 

It is interesting to compare the regions obtained by the above method with 
those obtained in Ref. [37], where the authors used the only viable mode for a low 
mass Higgs boson at the time {Mh = 120 GeV), 6677, to extract A from the visible 
mass distribution. After background subtraction, their best limit at 600 fb~^ was 
A e (0.26,1.94) at la. Here, for the 66r+r- we obtain A e (0.62,1.52), for the 
bbW+W- mode we obtain A e (0.41,2.20) and for the 6677 mode, A e (0.26,3.99). 
It is evident that the ratio provides a comparable exclusion region, especially consid- 
ering the fact that Ref. [37] considers relatively optimistic background subtraction. 
However, the ratio possesses advantages over the distribution analysis that may con- 
tain systematic uncertainties induced by the modelling of the shapes of both the 
signal and background. Note that an interesting study of the theoretical sensitivity 
of different initial states {gg — )■ HH, qq' — > HHqq', qq' WHH and qq ZHH) 
on the trilinear coupling can be found in [35]. 

It is important to consider the more general case when the top quark coupling to 
the Higgs boson is 7^ 1 [64-72]. This, for example, can be the case in a model where 
the top quark does not obtain its mass from a SM Higgs boson alone. However, one 
must keep in mind that additional heavy states potentially present in such models 
can further modify the predictions. A direct measurement of the magnitude of the 
Hti coupling can be made by observation of associated production of a single Higgs 
with top quark pairs [73] using boosted jet techniques that exploit the substructure 
of so-called 'fat' jets. Since the cross section for Higgs pair production, as well as the 
single Higgs cross section, both depend on the top coupling, yt, determination of yt 
and the triple coupling. A, cannot be done independently through measurement of the 
ratio Chh- Since the error on a determination of yt is expected to be 0{15%) [60], 
an investigation of the possible constraints in the yt — X plane is essential. This 
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Figure 6: The expected exclusion for A at one and two standard deviations for a given 
value of Atrue at 600 fb^^ for the bbT^T~ decay mode. The exclusion constructed from 
the ratio, Chh, is shown on the top panel, whereas the exclusion obtained from the cross 
section, anH, is shown on the bottom panel. 



can be done with a certain assumption on the underlying true value of A and i/t. 
For simplicity, here we assume that Atme = 1 and yt^truc = 1: i-e., they both have 
their respective SM values. We then calculate the induced error as we did before 
and calculate the la and 2a confidence levels on where the actual measurement will 
likely end up in the yt — X plane. The results are shown in Figs. 9, 10 and 11 for 
bbT~^T~ , bbW~^W~ and 6677 respectively, given an integrated luminosity of 600 fb~^. 
The figures illustrate an important point: for a model- independent determination of 
the Higgs triple self-coupling, a good measurement of yt is crucial. If, for example, 
we consider yt at the edges of the expected (9(15%) error, then yt = 0.85 yields 
A G (0.2, 1.03) whereas yt = 1.15 yields A G (1.18, ~ 2.5), both at la. This is a result 
of the sensitivity of the single and double cross sections on yt (given by Eq. (2.3)). 
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bb W yy- channel, y, =1, M„ =125 GeV, LHC@14 TeV, MSTW2008nlo68 cl 




-1.00 -0.75 -0.50 -0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 

^tme 



Figure 7: The expected exclusion for A at one and two standard deviations for a given 
value of Atrue at 600 fb~^ for the bbW~^W~ decay mode, constructed by using the ratio of 
cross sections (7//^^. 



bb77 channel, y, =1, M,^ =125 GeV, LHC(ai4 TeV, MSTVy2008nlo68cl 
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Figure 8: The expected exclusion for A at one standard deviations for a given value of 
Atrue at 600 fb^^ for the 6677 decay mode, constructed by using the ratio of cross sections 
Chh- The two standard deviations exclusion is not shown since it is weak. 

5. Conclusions 

We have considered the theoretical error on the ratio of cross sections of double-to- 
single Higgs production, Chhi at a 14 TeV LHC, including scale variation and par- 
ton density function uncertainties. Under the assumption that the double and single 
Higgs boson production cross sections possess a similar form of higher-order correc- 
tions, which we motivated in Section 3, we showed in the same section that the ratio 
is a more theoretically stable quantity than the cross section itself. Subsequently, 
assuming a 5% total theoretical error on Chh, and using conservative assumptions 
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bbr^T- channel, M„ =125 GeV, LHC@14 TeV, MSTW2008nlo68cl 
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Figure 9: The la and 2a confidence regions in the yt — ^ plane at 600 fb for the bbr^T 
decay mode, derived using Chh, under the assumption that Atrue = 1 and yt, true = 1- 
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Figure 10: The la and 2a confidence regions in the yt — ^ plane at 600 fb for the 
bbW~^W~ decay mode, derived using Chh, under the assumption that Atrue = 1 and 

yt, true = 1- 

on the experimental uncertainties of the quantities involved in measuring the ratio, 
we used this ratio to construct possible exclusions given a true value of the Higgs 
self-coupling at a 14 TeV LHC and integrated luminosities of 600 fb~^ and 3000 fb~^. 
In the SM (Atrue = 1 and yt,true = 1)? we concluded that it is possible to constrain 
the trilinear coupling to be positive, at 95% confidence at 600 fb~^, only using the 
discovery of the three viable channels. We also showed that a naive combination 
of the 'uncertainties' at Icr over the three channels indicates that a measurement of 
accuracy +30% and —20% is possible simply by using the ratio Chh at 3000 fb~^. 
The present work outlines the most precise method of determination of the Higgs 
triple self-coupling to date. We have also considered the more model-independent 
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bb77 channel, M„ =125 GeV, LHC@14TeV, MSTW2008nlo68cl 
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Figure 11: The la and 2a confidence regions in the yt — ^ plane at 600 fb for the 6677 
decay mode, under the assumption that Atrue = 1 and yt,tTue = 1- 

scenario where the top-Higgs coupling is not assumed to have the SM value, and 
we have considered the possible exclusion region in the yt — X plane in that case. 
Thus we concluded that an accurate, model-independent, determination of the Htt 
coupling, yt, is crucial to the determination of the Higgs boson triple self-coupling. 

We stress that the derived expected exclusion regions are based on the assump- 
tion that the SM is valid everywhere else except in the Higgs sector itself, and we 
assume no higher dimensional operators. Thus, deviations from expected exclusions 
would be an indication of some inconsistency in these assumptions that would require 
further assessment in the form of new physics models. Furthermore, it is obvious 
from the present study, as well as previous ones, that the measurement of the Higgs 
boson trilinear self-coupling is a challenging task, and further effort, both on behalf 
of theorists and experimentalists, should be made in order to obtain the best possible 
measurement during the lifetime of the LHC 
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